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Abstract 

The joint use of counting functions, Hilbert basis and Markov basis allows to define 
a procedure to generate all the fractions that satisfy a given set of constraints in 
terms of orthogonality. The general case of mixed level designs, without restrictions 
on the number of levels of each factor (like primes or power of primes) is studied. This 
new methodology has been experimented on some significant classes of fractional 
factorial designs, including mixed level orthogonal arrays. 
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1 Introduction 



All the fractional factorial designs that satisfy a set of conditions in terms 
of orthogonality between factors have been described as the zero-set of a 
system of polynomial equations in which the indeter minates are the com- 
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the reader. In Section [3] we write the problem of finding fractional factorial de- 
signs that satisfy a set of conditions as a system of linear equation s in which 



the in determinates are positive integers. In section HJ using 4ti2 (j4ti2 team 



( 120071 )) we find all the generators of some classes of fractional factorial de- 
signs, including mixed level orthogonal arrays and sudoku designs. Finally, in 
section [5] we consider the moves between different fractions as integer valued 
functions defined over the full factorial design. We build a procedure to move 
between fractions that use Markov basis. 



2 Notation and background 



2. 1 Full factorial design 



We adopt the notation used in iPistone and Rogantinl (120081 ) and denote: 



by Vj a factor with Uj levels coded with the n^-th roots of the unity: 

Vj = {ujQ, . . . , LJn^-i} Uk = exp (^i'^ kj ; 

by V the full factorial design with complex coding 

V = Vi X ■ ■ - Vj ■ ■ ■ X Vm . 

by ^T> the cardinality of V. 

by L the full factorial design with integer coding 

L = Z,„^ X ■ ■ ■ X Z„^, ■ ■ ■ X llrim 5 

by a an element of L 

a = (ai, . . .,am) = 0, . . . , - 1, j = 1, . . . , m . 
by [a — P] the m-tuple made by the componentwise difference 

(["1 - A]ni , • • • J«J - Pj]n, , • • • , [ttm - /?m]„„) ; 

the computation of the j-th element is in the ring Z„ . 
by Xj the j-th component function, which maps a point to its i-th compo- 
nent: 

Xj : P9 (Ci,...,Cm) ^ eP, ; 
the function Xj is called simple term or, by abuse of terminology, factor. 
by X" the interaction term X^^ ■ ■ -X^, i.e. the function 



X": v^{(^,...,Cm) ^ cr---C 



am . 

m 1 



We notice that L is both the full factorial design with integer coding and the 
exponent set of all the simple factors and interaction terms and a is both 
a treatment combination in the integer coding and a multi-exponent of an 
interaction term. 

The full factorial design in complex coding is identified as the zero-set in C"^ 
of the system of polynomial equations 

X;^-1 = , j = l,...,m. (1) 
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Definition 2.1 (i) A response / on the design V is a C-valued polynomial 
function defined on T). 

(2) The mean value onV of a response f , denoted by Exi{f), is: 

= ^ E /(C) • 

(3) A response f is centered on V if E-^lf) = 0. Two responses f and g are 
orthogonal on V if Ex>{f g) = 0, where g is the complex conjugate of g. 

It should be noticed that the set of all the responses is a complex Hilbert space 
with the Hermitian product: 

f-g = EM g) . 

Moreover 

(1) X"X^ = X["-'3]; 

(2) Ei,(XO) = 1, and ^©(X") = for a ^ 0. 

The set of functions {X" , a G L} is an orthonormal basis of the complex 
responses on design T>. In fact #L = ^T> and, from properties (i) and (ii) 
above, it follows that: 



EviX^'X'^) = Ev{X^"''^^) 



1 ifa = /3 
if a ^ /? 



In particular, each response / can be represented as a unique C-linear combina- 
tion of constant, simple and interaction terms. This representation is obtained 
by repeated applications of the re- writing rules derived from Equations ([1]). 
Such a polynomial is called the normal form of f on V. In this paper we intend 
that all the computation are made using the normal form. 

Example 2.1 Consider the 2^ full factorial design. All the monomial re- 
sponses on V are 

1, Xi, X2, X3, X1X2, X1X3, X2X3, X1X2X3 

or, equivalently, 

Xiom^ j^(i,o,o)^ j^(o,i,o)^ j^(o,o,i)^ j^(i,i,o)^ j^(i,o,i)^ X^^'^'^) 

and L is 

L = {(0, 0, 0), (1, 0, 0), (0, 1, 0), (0, 0, 1), (1, 1, 0), (1, 0, 1), (0, 1, 1), (1, 1,1)}. 



3 



2.2 Fractions of a full factorial design 

A fraction is a multiset (JF^,,/^) whose underlying set of elements JF^, is 
contained in V and /* is the multiplicity function /.j : JF^ — > N that for each 
element in JF^ gives the number of times it belongs to the multiset JF. 

All fractions can be obtained by adding polynomial equations, called gener- 
ating equations to the design equations [H in order to restrict the number of 
solutions. 

Definition 2.2 If f is a response on T> then its mean value on T, denoted 

by Ejr{f), IS 

EAf) = ^ E /(C) 

where #JF is the total number of treatment combinations of the fraction. 

A response f is centered if Ejr[f) = 0. Two responses f and g are orthogonal 
onTtfE^f g) = Q. 

With the complex coding the vector orthogonality of two interaction terms 
X" and as defined before (with respect to a given Hermitian product) 
corresponds to the combinatorial orthogonality (all the level combinations 
appear equally often in 

We consider the general case in which fractions can contain points that are 
replicated. 

Definition 2.3 The counting function R of a fraction T is a response defined 
on V so that for each ( R{C) equals the number of appearances of( in the 
fraction. A 0—1 valued counting function is called indicator function of a single 
replicate fraction T . We denote by Cq, the coefficients of the representation of 
RonV using the monomial basis {X", a E L}: 

i?(C) = Ec.^"(C) CeP gC. 

As the counting function is real valued, we have = C[_q,]. We will write cq 
in place of co,...,o- 

Remark 2.1 The counting function R coincides with multiplicity function 

Proposition 2.1 Let be a fraction of a full factorial design V and R = 
J2aeL Ca^" be its counting function. 
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(1) The coefficients Ca are: 



in particular, cq is the ratio between the number of points of the fraction 
and that of the design. 

(2) In a fraction without replications, the coefficients are related according 
to: 

/3eL 

(3) The term X" is centered on T , i.e. Ejc-(X°), if, and only if, 

Ca = C[_«] = . 

(4) The terms X"" and X^ are orthogonal on T , i.e. E_;f(X° X^) = 0, if, and 
only zf 

C[a-/3] = . 

Example 2.2 We consider the fraction T = {(—1, —1, 1), (—1, 1, —1)} of the 
2^ full factorial design of Example \2.1\ All the monomial responses on T and 
their values on the points are 



c 


1 




X2 


^3 


X1X2 


X1X3 


X2X3 


X1X2X3 


(-1,-1,1) 


1 


-1 


-1 


1 


1 


-1 


-1 


1 


(-1,1,-1) 


1 


-1 


1 


-1 


-1 


1 


-1 


1 



Using Item [I] of Proposition \2.1l it is easy to compute the coefficients Ca: 
C(o,i,o) = C(o,o,i) = C(i,i,o) = C(i,o,i) = 0; C(o,o,o) = C(i,i,i) = | and C(i,o,o) = 
C(o,i,i) = — |- Hence, the indicator function is 

F = - {1 — Xi — X2X3 + X1X2X3) . 

From the null coefficients we see that Xi and X3 are centered and that Xi is 
orthogonal to both X2 and X^ . □ 



2.3 Projectivity and orthogonal arrays 



Definition 2.4 A fraction T factorially projects onto the I-factors, I C 
{1, . . . , m\, if the projection is a multiple full factorial design, i.e. a full fac- 
torial design where each point appears equally often. A fraction T is a mixed 
orthogonal array of strength t if it factorially projects onto any I-factors with 
i^I = t. 
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Strength t means that, for any choice of t columns of the matrix design, all 
possible combinations of symbols appear equally often. 



Proposition 2.2 (Projectivity) (1) A fraction factorially projects onto 
the /-factors if, and only if, all the coefficients of the counting function 
involving only the I-factors are 0. 

(2) If there exists a subset J o/ {1, . . . , m} such that the J-factors appear in 
all the non null elements of the counting function, the fraction factorially 
projects onto the /-factors, with I = J'^ . 

(3) A fraction is an orthogonal array of strength t if, and only if, all the 
coefficients of the counting function up to the order t are zero: 

Cq = for all a of order up to t, a 7^ (0, 0, . . . , 0) . 

Example 2.3 (Orthogonal array) The fraction of a 2^ full factorial design 

.Fo = {(-1, -1, -1, -1, -1, 1), (-1, -1, -1, 1, 1, 1), (-1, -1, 1, -1, -1, -1), 

(-1,-1, 1,1, 1,-1), (-1,1, -1,-1, -1,-1), (-1,1, -1,1, 1,-1), (-1,1, 1,-1, 1,1), 
(-1,1, 1,1, -1,1), (1,-1, -1,-1, 1,1), (1,-1, -1,1, -1,1), (1,-1, 1,-1, 1,-1), 
(1,-1, 1,1, -1,-1), (1,1, -1,-1, 1,-1), (1,1, -1,1, -1,-1), (1,1, 1,-1, -1,1), 
(1,1,1,1,1,1)} 

is an orthogonal array of strength 2; in fact, its indicator function 



F — - + — -X1X4X5 + -XiX^X^Xq + -X1X3X4X5 

+ --^1-^2-^4-^5 + --^1-^3-^4-^5-^6 + --^1-^2-^4-^5-^6 

00 o 

+ g -^1-^2 -^3-^4-^5 ~ g -^1-^2-^3-^4-^5-^6 

contains only terms of order greater than 2, together with the constant term. 
□ 



3 Counting functions and strata 



From Proposition 12.11 and Proposition 12.21 we have that the problem of find- 
ing fractional factorial designs that satisfy a set of conditions in terms of 
orthogonality between factors can be written as a polynomial system in which 
the indeterminates are the complex coefficients of the counting polynomial 
fraction. 
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Example 3.1 Let's consider 3 factors, each one with two levels. The indicator 
functions F = J2a CaX^ such that the terms Xi, X2, X3 are centered on T and 
the terms Xi, Xj i,j = 1,2,3, i ^ j are orthogonal on T , where J-" = {( eV : 
F{() = 1}, are those for which the following conditions on the coefficients of 
F holds 

{Co = Cg + C123 
C123 = 2CoCi23 



Apart from the trivial F = 0, i.e. J-" = ^ and F = 1, i.e. J-" = V we find 
F = 1(1 + X1X2X3) anc? F = 1(1 - X1X2X3) 

Let's now introduce a different way to describe the full factorial design V and 
all its subsets. Let's consider the indicator functions 1^ of all the single points 
of P 

'1 C = (Ci, 

Cy^{Cl,...,U 

It follows that the counting function i? of a fraction can be written as 



Ic: V3iCi,...,Cm) ^ 



with = R{C) G {0,1, . . . ,n, . . .}. The particular case in which R is an 
indicator function corresponds to G {0, 1}. 

The coefficients are related to the coefficients Cq, as in the following Propo- 
sition 13.11 

Proposition 3.1 Let J-" be a fraction of V. Its counting fraction R can he 
expressed both as R = J2a (^aX" and R = ^^(^evyc^c- The relation between the 
coefficients Ca and yc^ is 

1 



Proof. From Proposition 12.11 we have 



□ 
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3.1 Strata 



As described in Section |2l we consider m fac tors, T>^ T>m wher e V 



Vtn = {^Oi ■ ■ ■ , -i}, for j = 1, . . . , m. From iPistone and RogantinI (120081 ) 



we recall two basic properties which hold true for the full design V 

Proposition 3.2 LetXj the simple term with level setQn = {^o, ■ ■ ■ -i}- 
Let's consider the term XJ and let's define 



1 r = 

nj/gcd{r, nj) r > 



Over V, the term Xj takes all the values of equally often. 

Proposition 3.3 Let X" = Xf ^ ■ ■ ■ X^™ an interaction. Xf ' takes values in 
where Si is determined according to the previous Proposition \3.^ Let's 
define s = lcm{si, . . . , s^). Over T>, the term X" takes all the values of VL^ 
equally often. 

Let's now define the strata that are associated to simple and interaction terms. 

Definition 3.1 Given a term X°, a G L = x . . . x the full design V 
is partitioned into the the following strata 

where uj^ G VL^ and s is determined according to the previous Propositions 1 3. 2\ 
and\33. 

Remark 3.1 We define strata using the conjugate X" of the term in place of 
the term X" itself because it will simplify the notations. 

Remark 3.2 Each stratum is a regular fr action whose defining equation is 



X"(^) = uj_h, Pistone and Rogantin (2008 ). 



We use na,h to denote the number of points of the fraction JF that are in the 



stratum D^, with h = 0, . . . , s — 1, 



na,h = yc 



The following Proposition 13.41 links the coefficients Ca with na,h- 

Proposition 3.4 Let J-' he a fraction ofV with counting fraction R = J2aeL Ca^"^ 
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Each G L, depends on na,h, h = 0, . . . , s — 1, as 



1 



where s is determined by X" (see Proposition \3. 3\) . Viceversa, each na^h,h = 
0, . . . , s — 1, depends on C[^ka], k = 0, . . . , s — 1 as 

na,h = X! '^[-ka]i^[hk] 

^ k=0 

Proof. Using Proposition 13.11 it follows that we can write the coefficients 
in the following way 

1 1 ■^-1 1 

TT fGf " h=G CeD'^ tF^ h=0 

For the viceversa, we observe the indicator function of strata can be obtained 
as follows. We define 

ifC^i 



fc=0 



s ifC = l 



We have FQ{uJk) = for all Uk & ^Ig, k ^ 0. It follows that 

Fa,o{o = -mn = -(i+c + ...+ c^^-'^" 

s s ^ 

is the indicator function associated to Dq. 

The indicator of = {( G P : X'^(C) = uJh} = {( eV : X"(C) = ui-h]} will 
be 

We get 

n„,,= E i?(C) = 5: F,,,(C)/?(C) = 

= 2^ ^[kh]Cp = 2^ i^[kh]C[-ka] 

^ k,l3:[ka+f3]=0 ^ k=0 



□ 
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Remark 3.3 From Proposition\3.4\ we get 



no,h = 0, h = 1, . . . ,s -1 

na,0 = 2^ C[_fcQ,] 

^ k=0 

and in particular uq^q = #JF. 



We now use a part of Proposition 3 of iPistone and RogantinI (120081 ) to get 
conditions on /j that makes centered on the fraction JF. 

Proposition 3.5 Let X"" he a term with level set Vtg on full design T>. Let 
P{C) the complex polynomial associated to the sequence {na,h)h=o,...,s-i so that 

h=0 

and let's denote by $s the cyclotomic polynomial of the s-roots of the unity. 

(1) Let s he prime. The term X" is centered on the fraction T if, and only 
if, its s levels appear equally often: 

nafi = n-a,! = . . . = ria^s-i = 

(2) Let s = pi^ . . .p^'' withpi prime, fori = 1, . . . ,d. The term X" is centered 
on the fraction T if, and only if, the remainder 

H{0 = P{Q mod $,(C) 

whose coefficients are integer linear combinations ofua^h, h = 0, . . . , s — 1, 
is identically zero. 



Proof. See Proposition 3 of iPistone and RogantinI (120081 ). □ 

Remark 3.4 Being Df[ a partition ofV, if s is prime we get Aq, = 

If we remind that n^^h are related to the values of the counting function R of 
a fraction JF by the following relation 

na,h = Vo 

this Proposition 13.51 allows to express the condition X" is centered on T as 
integer linear combinations of the values -R(C) of the counting function over 
the full design D. In the Section HJ we will show the use of this property to 
generate fractional factorial designs. 
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We conclude this section limiting to the particular case where all factors have 
the same number of levels s and s is prime. We provide some results concerning 
the coefficients of counting functions, regular fractions, wordlength patterns 
and margins. 



3.2 Coefficients of the polynomial counting function 



From Proposition 13.51 we get the following result on the coefficients of a count- 
ing function 

Proposition 3.6 Given a counting function R = I]q,Cq,X", if Ca = then 
C[k-a] = for all k = 1, . . . , s — 1, where [k ■ a] is a + . • . + a in the ring Z™. 

k times 

Proof. Let's consider Ck-a- From Proposition 13. 5[ Ck-a is equal to zero if, and 
only if, 

E yc= yc = ■■■= E 

We observe that 



where [kh] is h + . + h in the ring Z^. 

k times 

It follows that X" and X^'" partition V in the same strata and therefore we 
get the proof. □ 



3.3 Regular designs 



Let's consider a fraction wit hout replicates and wi t h ind icator function 
F = Z^aCa^"- Proposition 5 in (jPistone and Rogantinl (120081 ) ) states that a 
fraction JF is regular if, and only if, its indicator function F has the form 



where £ C L, £ is a subgroup of L and e : C 
mapping. 



{uq, . . .,ujs-i} is a given 
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If we use Proposition 13.51 we immediately get a characterisation of regular 
fractions based on the frequencies Ua^h- 

Proposition 3.7 Given a single replicate fraction T with indicator function 
F = J2a (^aX^ the following statements are equivalent: 

(i) T is regular 

(a) forria^h there are only two possibilities 

(a) if Ca = then na,h = h = 0, . . . , s — 1, 

(b) if Ca 7^ then 3/i* G {0, . . . , s — 1} such that 



na,h 



^ ifh = K 
otherwise 



Proof. Using Proposition 13.41 we get 



1 



Proposition 5 in iPistone and Rogantinl (120081 ) gives the following conditions 
on the coefficients of the indicator function F of a regular fraction JF: 



e(") 
I • 

otherwise 



aeCCL 



where e : C {uq, . . . , Us^i}, I = jj^mathcalL and £ is a subgroup of L. 
Let's consider a E C. We get 



1 e{a) 

Let's suppose e(a) = ujh,. We obtain 

na,hUJh + {^na,h. - hujh, = (2) 
To simplify the notation we let ah = -^ria^hi h = 0, . . . , s — 1, h h^, and 



ah. 



hj^a^h^ — J. Therefore, from the proof of item (1) of Proposition 13. 5[ 



for the relation [2] to be valid, it should be 

aQ = ai = . . . = tts-i 
Being Y^I^Jq na,h = #^ it follows 

E^«,/.= E (#^^)a, + (#I))(a,. + y) = (#I?)X:a. + ^ = #^ 

h=0 h=Q,h^K ' h=0 ' 
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and so 

We finally get 



na,h 



otherwise 



Being C a subgroup of L it follows that G £ and so cq = 1//. We also know 
that Co = ^ and therefore 

For the null coefficients of F, {cq, : a G L — £}, it is enough to use Proposition 
?? to conclude the proof. □ 



3.4 Wordlength Pattern 



Aberration is often used as a criterion to compare fra ctional facto r ial de - 
signs. The generalized minimum aberration, proposed bvlXu and Wul f|2001l). 
i s base d on the generalised wordlength pattern, see also lBeder and Willenbring 
(120091 ). It can be shown that the generalized wordlengths can be written in 
terms of the squares of the modules of the coefficients c^, obtaining 



A, 




for j 



, m 



''0 wt{a)=j 



where wt{a) is the Hamming weight of a, i.e. the number of nonzero compo- 
nents of a. We now express the square of the module of the coefficient Cq, in 
terms of Ua^h- 

Proposition 3.8 



Y.i'^Xh - na,hn[a,h~^]) /or 7 G {1, . . . , s - 1} 

h=0 



Proof. From Proposition 13.41 we get 



1 '-^ 



It follows 
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s-1 s-1 



Krt^) h=Q k=o 



s—l s—1 



vfr'^) h=Q k=0 



s—l s—l 



7=0 p=0 



2 

\ca\ must be a real number. Being = 1 it follows 



2 s—l 2 
( (^pj2 ^p^^'P ~ |Ca|^)^0 + ^^pp E E ^".p'^I^.P-tFt = (3) 

To simplify the notation we let Oq = ( (^p)2 I]p=o ^a,p " = 
1^)2 I]p=o ^a,p^[a,p-7] ) 7 = 1, . . . , s — 1. Therefore, by Lemma ??, for the re- 
lation [3] to be valid, it should be 

= Ctl = . . . = Cls-l 

Using one of the equalities, ao = ah h = 1, . . . , s — 1, it follows 



s-1 

C, 



1' = ^ 



\2 



□ 



Remark 3.5 Proposition \ 3.S\ provides a useful tool to compute the modules 
of the coefficients Ca- Indeed it is enough to choose 7 = 1 and compute \caf' 



as 



\))2 Z]h=o(^a,h f^ot,hfT'[oL,h-l]) 



Remark 3.6 We make explicit these relations for 2 and 3 level fraction. 
If s = 2 then 



I |2 l / \2 



If s = 3 then, choosing 7 = 1, 

I |2 _ ^ i 2,2,2 \ 

Remark 3.7 We observe that, denoting byfia the mean of the values ofna,h, 
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na = lT,lJona,h, we get 



s-l s-1 

,2 



h=0 h=0 

We have 



1 



n 



■2 

h,k=0 



2 ^ ^ ^a,h'^a,k 



/s-l s-l s-l \ 

^ \h=0 h=0 h=0 J 

where s^, = Proposition \3.8\ states that all the quantities Z]h=o ^o,h^a,[/i-7] 
are equal and so, choosing, without loss of generality, j = 1, we get 

/s-l s-l \ /s-l s-l 

"^1 = — (Yl ^l,h + 2s* Y ^a,hna,lh-i] ) = — (Y ^l,h + (s - 1) II n^,hna,[h-i] 

^ \h=0 h=0 / ^ \h=0 h=0 

and therefore 



s-l s-l 

,2 —2 



Y (^a./^ "~ ^")^ = Y 



h=0 h=0 

s - 1 /'-'^ \ 
= —— Y ^l,h - Y ^o.,hna,yh-i] = 

^ \h=0 h=0 ) 

s 

It follows that, if we denote by a"^ the variance ofua^h, = j Z]h=o ~ ^«)^ 
we get 



i2 



2 



and so the square of the module of Cq represents, apart from a multiplicative 
constant, the variance of the frequencies ria^h- 



3. 5 Margins 



We now examine the relationship bet ween the margins and the co efficients 
of the counting functions. We refer to (jPistone and Rogantin and we 

report here a part of it. 
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For each point ^ G P we consider the decomposition ( = {(j, (j) where I C 
{1, . . . , m} and J = {1, . . . , m} — I = P is its complement. We denote by 
RiiCi) the number of points in JF whose projection on the I factors is (j. 

In particular if J = {1, ... , m} we have Rj = R and if / = we have Rj = 

We denote by Lj the subset of the exponents restricted to the I factors and 
by a I an element of Lj: 

Li = {aj = (ai, . . . , a^), aj = if j e J} 

Then for each a E L and E V we have a = aj + aj and = 
Xf{(j)Xj'{(j). Finally we denote by Vj and T>j the full factorial over the 
/ factors and J factors, respectively {V = Vj x T>j). 

We have the foUowirig prop osition (see item 1 and 2 of Proposition 4 of 



Pistone and RogantinI (120081 ) ) 



Proposition 3.9 Given a fraction ofD 

(1) the number of replicates of the points of T projected on the I factors is: 

RiiCi) = #^jT.^»,X'^'iCi) 

(2) T fully projects on the I factors if, and only if, 



Rj{Cj) = #Vj ■ Co = 



We will refer to Rj as /c-margin, where k = #J. The number of /c-margins is 
(j^^ and each fc-margin can be computed over points (j G Vj. It follows 
that there are (1 + s)™ marginal values in total. 

Using item 1 of Proposition 13.91 and reminding that we work with a prime 
number of level s we have 



or, by the definition of Rj as the restriction of R over the / factors. 



We point out the following relationship between margins. 

Proposition 3.10 If A C B C {!,..., m} and Rb{Cb) = s'^'^'^cq then 
Ra(Ca) = s'^-^^Cq where #S = A;^ and #A = kA 
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Proof. Let's put Ai = B — A. We have 

RAiCA)= E RauaACa,CaJ= E RBiCAXAj=s'^-'^S^-'-Co = s'^-'-Co 
□ 

We finally observe that, as we already pointed out, given C C L a set of 
conditions = 0, a G C translates in a set of conditions SceD^ Vc = h = 
0, . . . , s — l,a G C where A does not depend by a (and by h). In general, 
with respect to margins, the situation is different. For example let's suppose 
to have a JF that fully projects over the Ji and the I2 factors, with Ji fl /2 = 
and #/i 7^ 7^/2. From Proposition 13.91 we obtain 

RhiCh) = ^77 and Ri^iCh) = ^TfJ 



4 Generation of fractions 



Let use strata to generate fractions that satisfy a given set of constrains on 
the coefficients of their counting functions. Formally we give the following 
definition 

Definition 4.1 A counting function R = J2aCaX°' associated to J-" is a C- 
compatible counting function if its coefficients satisfy to 

c„ = 0, a G C, C C Z„j X . . . 

We will denote by OF{ni . . . rim, C) the set of all the fractions whose counting 
functions are C-compatible. 

In the next sections, we will show our methodology on Orthogonal Arrays and 
Sudoku designs. 



4.1 OA{n,s"',t) 



Let's consider OA{n, s^, t), i.e. orthogonal arrays with n rows and m columns 
where each columns has s symbols, s prime and with strength t. 

Using Proposition 12 . 21 we have that the coefficients of the corresponding count- 
ing functions must satisfy the conditions Cq = for all a G C where C C L = 
{a : < ||a;|| < t} where ||a|| is the number of non null elements of a. We 
have Ni = J2k=i [kji-^ " 1)^ coefficients that must be null. 
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It follows that OF(s™,C) = [jnOA{n, s"',t). 

Now using Proposition 13.51 we can express these conditions using strata. If we 
consider a G C we write the condition = as 



To obtain all the conditions it is enough to vary a E C. We use Proposition [XB] 
to limit to the a that give different strata. It is easy to show that we obtain 
N2 = different a, each of them generate s linear equations, for a total of 

constraints on the values of the counting function over P. 
We therefore get the following system of linear equations 



AY = XI 



where A is the (A^ x s™') matrix whose rows contains the values, over T>, of the 
indicator function of the strata, 1 no , Y is the column vector whose entries 

are the values of the counting function over T), X will be equal to ^ and 
1 is the s*" column vector whose entries are all equal to 1. We can write an 
equivalent homogeneous system if we consider A as a new variable. We obtain 

Ay = 



where 



A 



A 



[A -1] 



and 



Y 



Y 



iY,X) 



In an equivalent way, we can also express the conditions Cq = for all a G C 
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in terms of margins. We obtain 

RiiCi) = s^'-^*'^co 

where / C {1, . . . , m} and 1 < 7^/ < t. If we recall Proposition I3.10[ we can 
limit to the margins Rj where 7^/ = t. We have values of such t margin 

In this case, with the same approach that we adopted for strata, we obtain a 
system of linear equations 

BY = pi 

where p = s^~^Cq and its equivalent homogeneous system 

BY = 

Now we can find all the generators of OF{s"^,C, that means of Orthogonal 
Arrays OA{n, s"^,t), by computing the Hilbert B asis corresponding to A (or , 



equivalently, to B). This approach is the same of ICarlini and Pistond (120071 ) 
but, in that work, the following conditions were used 



The advantage of (or margins) is that we avoid computations with 

complex numbers (X"(^)). W e explain this po int in a couple of examples. For 
the computation we use 4ti2 fl4ti2 team! fl2007h ). 



We use both A (strata) and B (margins) because, even if they are fully equiv- 
alent from the point of view of the solutions that they generate, they perform 
differently from the point of view of the computational speed. 



4.1.1 OA(r^,2^2) 



0A{n,2^ ,2) were investigated in ICarlini and Pistond (120071 ). We build both 
the matrix A and B. They have 30 rows and 40 rows, respectively and 33 
columns. We find the same 26, 142 solutions as in the cited paper. 



4.1.2 OA{n,3^,2) 

We build both the matrix A and B. They have 54 rows and 27 rows, respec- 
tively and 28 columns. We find 66 solutions, 12 have 9 points, all different and 
54 have 18 points, 17 different. 
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Finally we point out that 4ti2 allows to specify upper bounds for variables. 
For example, if we use B and we are interested in single replicate orthogonal 
arrays, we can set 1 as the upper bound for |/^, ( E V. The upper bound for 
the variable p can be set to s™~* = 3^-^ that corresponds to Cq = 1, i.e. to the 
full design V. 



4-2 OA{n,ni . . .nm,t) 



Let's now consider the general case in which we do not put restrictions on the 
number of levels. 



4.2.1 OA(n,42,l) 

In this case the number of levels is a power of a prime, 2^. Using Proposition 
12. 21 we have that the coefficients of the corresponding counting functions must 
satisfy the conditions = for all a G C where C C L = {a : = 1}. 

Let's consider ci q. From Proposition 13.21 we have that Xi takes the values in 
Qg where s = 4. From Proposition 13. 5^ Xi will be centered on JF if, and only 
if, the remainder 

HiC) = PiC) mod $4(0 
is identically zero. We have $4(C) = 1 + (see iLand (119651 )) and so we can 
compute the remainder 

^(C) = "'(1,0),0 - "'(l,0),2 + ('^(l.O),! - '^(1,0),3)C 

The condition H{() identically zero translates into 



^(l,0),0 - '^(l,0),2 = 
'^(1,0), 1 — "-(1,0), 3 = 

Let's now consider C2,o- From Proposition 13 . 21 we have that Xf takes the values 
in Qg where s = 2. From Proposition 13.51 Xf will be centered on JF if, and 
only if, the remainder 

HiC) = P(C) mod $2(0 
is identically zero. We have $2(0 = 1 + C (see Lang ( 19651 )) and so we can 
compute the remainder 



^(0 



n 



(2,0),0 — "-(2,0), 1 



If we repeat the same procedure for all the a such that ||a;| 
that 

na,h = J2 yc 



1 and we recall 
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orthogonal arrays Oy4(n, 4^,1) become the integer solutions of the following 
integer linear homogeneous system 



2/00 
2/10 

2/20 
2/30 

2/01 
2/11 

2/21 
2/31 
2/02 
2/12 
2/22 
2/32 
2/03 
2/13 
2/23 
2/33 



Using 4ti2 we find 24 solutions that correspond to all the Latin Hypercupe 
Designs (LHD). 
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4.2.2 0A{n,6'^,l) 

As in the previous examples, using Proposition l2.2l we have that the coefficients 
of the corresponding counting functions must satisfy the conditions = for 
all a G C where C C L = {a : ||a|| = 1}. 

Let's consider ci^o- From Proposition 13.21 we have that Xi takes the values in 
fig where s = 6. From Proposition 13.51 Xi will be centered on JF if, and only 
if, the remainder 

HiC) = PiC) mod $6(0 
is identically zero. We have $6(C) = 1 - C + (see IhaneJ fll965h ) and so we 
can compute the remainder 



H{C) = «(1,0),0 - '^(l,0),2 - '^(1,0),3 + «(1,0),6 + ('^(1,0),1 + '^(1,0),2 " '^(1,0),5 - ^(1,0),6)C 
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If we repeat the same procedure for all the a such that ||a|| = 1 and we recall 
that 

na,h = J2 

orthogonal arrays OA{n, 6^, 1) become the integer solutions of an integer linear 
homogeneous system AR = where the matrix A is built as in the previous 
case of 0A(n,4^, 1). Using 4ti2 we find 620 solutions that correspond to all 
the Latin Hypercupe Designs (LHD). 



4-3 Sudoku designs 



As shown in iFontana and Rogantinl (120081 ). a sudoku can be described using 



its indicator function. Here we report a very short synthesis of Section 1.3 of 
that work. 

A X with p prime sudoku design can be seen as a fraction of the full 
factorial design T): 

V = RiX R2XC1XC2X SiX S2 

where each factor is coded with the p-th roots of the unity. Ri and R2, Ci and 
C2, Si and S2, represent the rows, the columns and the symbols of the sudoku 
grid, respectively. 



The following proposition (Proposition 5 of Fontana and Rogantin ( 20081 )) 
holds. 

Proposition 4.1 Let F be the indicator function of a fraction J-' of a design 
design, F = J^a^L^a^'^ ■ The fraction T corresponds to a sudoku grid if and 
only if the coefficients satisfy the following conditions: 

(1) 6000000 = i-^- the ratio between the number of points of the fraction 
and the number of points of the full factorial design is l/p^; 

(2) for allij G {0, 1, . . . ,p - 1}.- 

(a) bi^i^i^i^oo = for (ii, ^2, «3, ^4) 7^ (0,0,0,0), 

(b) h^i^ooirje = for {ii,i2, is, ie) (0, 0, 0, 0), 

(c) booi^i^irje = for (i3,i4,«5,«6) 7^ (0,0,0,0), 

(d) 6»i0i30»5i6 = for iii.i^.i^.io) ^ (0,0,0,0) 

i.e. the fraction factorially projects onto the first four factors and onto 
both symbol factors and row /column/box factors, respectively. 

From this Proposition, we define C as the union of Ci, C2, C3 and C4, where 
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Ci = {(?iW400) : (21,22,^3,^4) ^ (0,0,0,0)} 
C2 = {iiii200i5i6) : (21,22,^5,%) (0,0,0,0)} 
C3 = {(OOisui^ie) : (zg, ^4, 25, ig) 7^ (0,0,0,0)} 
C^ = {(i^0i^0i^iQ) : (i;L,i3,i5,iig) ^ (0,0,0,0)} 



The problem of finding Sudoku becomes equivalent to find C-compatible count- 
ing functions, that are (i) indicator functions and (ii) that satisfy the additional 
requirement = l/p^. 



4.3.1 4x4 Sudoku 

We use the conditions C to build both the matrices A and B. A has 78 rows. 
With respect to B, that corresponds to the margins that must be constant, 
if we recall Proposition 13.101 we obtain 64 constraints, all corresponding to 
4-margins. 

To find all sudoku we use 4ti2, specifying the upper bounds for all the 65 
variables. The upper bounds for y^, ^ G must be equal to 1. If we use A, 
the upper bound for A must be set equal to ^ = y = 8, while if we use b 
the upper bound for p must be set equal to s™~'^6o = = 1. 



We find all the 288 different 4x4 sudoku as in lFontana and RogantinI ( 120081 ). 
We point out that to solve the problem using A the total time was 31.59 
minutes, while using B the total time was only 58.04 seconds on the same 
computer. 

If we admit counting functions with values in {0, 1, 2} and t^JF < 32 we find 
55, 992 solutions. 



5 Moves 



Sometimes, given a set of conditions C we are interested in picking up a solution 
more than in finding all the generators. The basic idea is to generate somehow 
a starting solution and then to randomly walk in the set of all the solutions 
for a certain number of steps, taking the arrival point as a new but still C- 
compatible counting function. 

Let's use the previous results on strata to get a suitable set of moves. We will 
show this procedure in the case in which all the factors have the same number 
of levels s, S prime, but it can also be applied to the general case. In Section H] 
we have shown that counting functions must satisfy the following set of linear 
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equations 



AY = XI 



where A corresponds to the set of conditions C written in terms of strata. 

It follows that if, given a C-compatible solution Y, such that AY = XI, we 
search for an additive move X such that A(Y + X) is still equal to Al, we 
have to solve the following linear homogenous system 

AX = 

with X = (a;^), C G P, G Z and + > for all C We observe that 
this set of conditions allows to determine new C-compatible solutions that give 
the same X. We know that A = — so this homogenous system determines 
moves that do not change the dimension of the solutions. 

Let's now consider the extended homogeneous system, where A has already 
been defined in Section HI 

Ax = 

with X = C ^ 1^, e Z and + X(; > for all ( &T>. 

Given Y = (y. Ay), where Y is C-compatible counting function and Ay = 

— ^ — , the solutions of AX = determine all the other Y + X = {Y + X, Xy+x) 
such that A{Y + X) =0. Y + X are C-compatible counting functions whose 
sizes, sXy+x, are, in general, different from that ofY. 



5.1 Markov Basis 



We use the theory of Markov basis (see for example iDrton et al.l (120091 ) where 
it is also available a rich bibliography on this subject) to determine a set of 
generators of the moves. 

We use the following procedure in order to randomly select a C-c ompatible 



counti ng function. We compute a Markov basis of ker(y4) using 4ti2 (j4ti2 team 



( 120071 )). Once we have determined the Markov basis of ker(y4), we make a 
random walk on the fiber of F, where Y , as usual, contains the values of 
the counting function of an initial design JF. The fiber is made by all the C- 
compatible counting functions that have the same size of JF. The randow walk 
is done randomly choosing one move among the feasible ones, i.e. among the 
moves for which we do not get negative values for the new counting function. 

In the next paragraphs we consider moves for the cases that we have already 
studied in Section HI 



24 



5.2 Orthogonal arrays 



5.2.1 Ov4(r^,2^2) 

We use the matrix A, already built in Section 14.1.11 and give it as input to 
4ti2 to obtain the Markov Basis, that we denote by Ai. It contains 5.538 
different moves. Given M = (x/^) & A4 we define = max(a;^, 0) and 
M- = max(-x^, 0). We have M = M+ - M". 

As an initial fraction JFq, we consider the eight-run regular fraction whose 
indicator function Rq is 

Ro = \{l + XiX2X3)(l + X1X4X5) 

We obtain the set of feasible moves observing that a move M G to be 
feasible, should be not negative when Ro is equal to zero that means 

(1 - Ro)M- = 

We find 12 moves. Analogously an element M ^ Ai such that 

(1 - Ro)M+ = 

gives a feasible move, —M. In this case we do not find any of such element. 

Therefore, given Rq, the set of feasible moves becomes J^Rq that contains 
12 + different moves. 

We randomly choose one move M^^ out of the 12 available ones and move to 

Ri = Ro + Mr, 

We run 1.000 simulations repeating the same loop, generating Ri as Ri = 
Ri_i + MR^_^. 

We o btain all the 60 diffe r ent 8- run fractions, each one with 8 different points 



as m 



Carlini and Pistond (120071 ). 



Using A we obtain the set M. that contains 18 different moves. 



5.2.2 OA(n,33,2) 

Using A as built in the Section I4.1.2[ we use 4ti2 to generate the Markov 
basis corresponding to the homogeneous system AX = 0. We obtain Ai that 
contains 81 different moves. 
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As an initial fraction we can consider the nine-run regular fraction whose 
indicator function Rq is 

Rq = -{1 + X1X2X3 + X^XgXg) 

We run 1.000 simulations repeating the same loop, i.e. generating Ri as Ri = 
R,-i + Mr^_^. 

We obtain all the 12 different 9-run fractions, each one with 9 different points 
as known in the literature and as found in Section I4.1.2[ 

Using A we also obtain the set M. that contains 10 different moves. 
5.2.3 4x4 sudoku 

Using the matrix A built in Section 14.3.11 we run 4ti2 getting the Markov 
basis M. that contains 34.920 moves. 

We randomly choose an initial sudoku 



3 2 

4 1 


4 1 

3 2 


2 3 
1 4 


1 4 

2 3 



The corresponding indicator function is 

Fo = ^(1 - i?2Ci5i^2)(l - R1C2S1) . 

Then we extract from M. the feasible moves. We obtain a subset M.Fo that 
contains 5 different moves. We repeat the procedure on — 7V1 and we obtain 
other 9 moves. 

We randomly choose one move Mp^ out of the 5 + 9 available ones and move 
to 

We run 1.000 simulations repeating the same loop Fi = Fj_i + Mi?._j. 
We obtained all the 288 different 4x4 sudoku. 
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6 Conclusions 



We considered mixed level fractional factorial designs. Given the counting 
function of a fraction T we translated the constraint Cq = 0, where Cq, is a 
generic coefficient of its polynomial representation R = Cq,X°, into a set of 
linear constraints with integer coefficients on the values that R takes on all 
the points ( E V. We obtained the set of generators of the solutions of some 
problems using Hilbert Basis. We also studied the moves between fractions. 
We characterized these moves as the solution of a homogeneous linear system. 
We defined a procedure to randomly walk among the solutions that is based on 
the Markov basis of this system. We showe d the procedure o n some examples. 
Computations have been made using 4ti2 fl4ti2 team! tOO^ ). 



Main advantages of the procedure are that we do not put restrictions on the 
number of levels of factors and that it is not necessary to use software that 
deals with complex polynomials. 

One limit is in the high computational effort that is required. In particular 
only a small part of the Markov basis is used because of the requirement that 
counting functions can only take values greater than or equal to zero. The 
possibility to generate only the moves that are feasible could make the entire 
process more efficient and is part of current research. 
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